Exact Results for Average Cluster Numbers in Bond Percolation on Lattice Strips 
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We present exact calculations of the average number of connected clusters per site, {k), as a 
function of bond occupation probability p, for the bond percolation problem on infinite-length 
strips of finite width L y , of the square, triangular, honeycomb, and kagome lattices A with various 
boundary conditions. These are used to study the approach of (k), for a given p and A, to its value 
on the two-dimensional lattice as the strip width increases. We investigate the singularities of (k) in 
the complex p plane and their influence on the radii of convergence of the Taylor series expansions 
of (k) about p = and p = 1. 



I. INTRODUCTION 



The study of percolation gives insight into a number 
of important phenomena such as the passage of fluids 
through porous media and the effect of lattice defects 
and disorder on critical phenomena. Here we consider 
bond percolation. Let G — G(V, E) be a connected 
graph defined by a set V of vertices (sites) and a set 
E of edges (bonds) connecting pairs of vertices. Multiple 
bonds connecting the same pair of vertices are allowed, 
although most of the graphs that we study have only sim- 
ple bonds. We denote the number of vertices and bonds 
as n = n(G) = \V\ and e(G) = \E\. One now envisions 
generating a large ensemble of corresponding graphs in 
which each bond is randomly present with probability 
p£ [0,1]. In the usual statistical mechanics context, one 
is interested in a <i-dimensional thermodynamic limit of 
a regular lattice graph A in which the bonds are present 
with probability p. Consider the probability P(A,p) that 
a given site belongs to an infinite cluster. For a given 
lattice A, as p decreases from 1, P(A,p) decreases mono- 
tonically until, at a critical value, p c ,A, it vanishes and 
remains identically zero for < p < p c .A- This has an 
important effect on physical phenomena that take place 
on such a bond-diluted lattice. For example, consider a 
spin-model (with finite-range spin-spin interactions) on 
A and assume that for the undiluted lattice the model 
is above its lower critical dimensionality so that there 
is a phase transition at a finite temperature T c . Then 
as p decreases below unity, T c (p) decreases below T c (l), 
and as p decreases to p c ,A> T c {p) vanishes nonanalytically 
and remains zero for < p < p c .A- Other quantities also 
behave nonanalytically at p c ,A, such as the mean size 



S(A,p) of the percolation cluster, which diverges as p in- 
creases through p c ,A- Some previous literature on bond 

fiercolation relevant to our present work includes Refs. 

An interesting quantity in this context is the number of 
connected components (clusters), including single sites, 
for a given lattice A, divided by the number of sites on the 
lattice and averaged over all of the graphs in the above 
ensemble. We denote this mean cluster number per site 
as (fc)A- Explicit exact values for (k)\ at the respective 
values p = p c ,A, denoted (k) c ,A, have been calculated 
in Ref. [22j for the square, triangular, and honeycomb 
lattices (where they were denoted nf~ A ). However, to 
our knowledge, (k) has never been calculated exactly as 
a function of p for the full range p £ [0, 1] for any lattice 
except a one-dimensional or Bethe lattice 0, El • 

In this paper we present exact calculations of this aver- 
age cluster number per site, (k), as a function of p, for a 
variety of infinite-length, finite- width strips of regular lat- 
tices. We consider strips of the square, triangular, honey- 
comb, and kagome lattices. These are of interest since, at 
least for modest strip widths, one can obtain explicit an- 
alytic expressions for (k) and can exactly determine, e.g., 
singularities that these expressions have in the complex p 
plane and their influence on series expansions. When re- 
ferring to a specific lattice strip A s (including transverse 
boundary conditions), we shall denote the average cluster 
number per site as (fe)A„ • Our results interpolate between 
the known exact solutions for the one-dimensional lattice 
(line) and the case of two dimensions, and complement 
numerical simulations and series expansions. Early stud- 
ies of percolation on infinite-length, finite-width strips 
include Refs. 0,0], which focused on universal proper- 
ties such as critical exponents. While the average cluster 
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number is obviously nonuniversal, depending as it does 
on the specific type of lattice, it still give useful informa- 
tion about the phenomenon of percolation, as we shall 
show. 

We take the longitudinal and transverse directions to 
be x and y and denote the size of the lattice strips in 
these directions as L x and L y and the respective bound- 
ary conditions as BC X and BC y . We focus on the limit of 
infinite length, L x — > oo, for which the results are inde- 
pendent of the longitudinal boundary conditions. For an 
infinite-length strip of a lattice A, as the width L y — > oo, 
one expects (fc) to approach a limiting function of p 
which is independent of the transverse boundary con- 
ditions and is equal to (k) for the corresponding infi- 
nite two-dimensional lattice A. In particular, for a given 
infinite-length, finite-width strip of the lattice A, it is 
of interest to evaluate our exact expressions for (k) at 
V = PcA and study how the resultant value approaches 
the critical value (fc) c ,A for the corresponding infinite two- 
dimensional lattice. 



II. CALCULATIONAL METHOD 

For a given graph G — (V, E) we calculate (A;) by mak- 
ing use of the fact that it can be expressed as a cer- 
tain derivative of the (reduced) free energy of the Potts 
model. We first recall this expression. Define a spanning 
subgraph G' = (V, E') as a subgraph of G with the same 
vertex set V and a subset of the edge (bond) set, E' C E. 
From the definition of the average number of clusters per 
site, we have 



(k) 



(!/") Sc k(G')p e(G ">(l - P )< G )-< G ') 
Ec P e{G ' ] (! ~ p)< G )-< G ') 



(l/n)E G ,fc(G>/(l--p)] e(G ' ) 
£ G '[p/(l-p)] e(G,) 



(2.1) 



This follows because each G 1 contains k(G') connected 
components, and appears in the numerator of the ex- 
pression in the first line with weight given by p e ^ G )(1 — 
p)e(G)-e(G )^ smce ^fie probability that all of the bonds 
in G' are present is p e ^ G ' and the probability that all 
of the other e(G) — e(G') bonds in G are absent is 
(1 — j>) e ( G )~ e ( G ). This sum in the numerator over the set 
of spanning subgraphs G' is normalized by the indicated 
denominator and by an overall factor of 1/n to obtain 
the average number of connected components (clusters) 
per site. As noted, we shall focus on the limit in which 
the strip length L x — > oo and hence n — > oo. 

The cluster representation for the partition function of 



the g-state Potts model is 0, Q 

Z(G,q,v)= J2 <Z fc(G V< G '> 



(2.2) 



G'CG 



In the Potts spin model, the spin Oi at each vertex i can 
take on q different values cr, = 1, 2, .., q and 



v = e K - 1 , K = $ J 



(2.3) 



where (3 — l/(fcsT) with T the temperature and J the 
spin-spin exchange constant. The formula 12.21 allows 
one to generalize q from the positive integers to a non- 
negative real variable while retaining a Gibbs measure 
for Z(G,q,v). 

On a finite graph G one defines the (reduced) free en- 
ergy per site of the Potts model as 



f(G,q,v)=HZ(G,q,v) 1 / n ] 



and, in the limit 



f({G},q,v)= lim f(G,q,v) 



(2.4) 



(2.5) 



where we use the symbol {G} to denote the formal limit 
linin^oo G for a given family of graphs. (The actual free 
energy per site is F = —ksTf; henceforth we shall simply 
refer to / as the free energy.) If one sets 



P 



(2.6) 



differentiates f(G, q, v p ) with respect to q, and sets q = 1, 
one obtains precisely the expression for (k) n , as given in 
eq. (j2.1jl . i.e., 



df(G, q,v p ) 



(*>» = 



In particular, in the limit n — > oo, 



9=1 



df({G},q,v p ) 



(2.7) 



(2.8) 



(k) 



Our method for calculating (fc) is to use eq. (I2.8|l in 
conjunction with exact results that we have computed 
for the free energy of the Potts model on infinite-length, 
finite- width strips of various lattices One can 

also study quantities in the percolation problem in an 
analytic manner without making use of eq. I|2.8(l (e.g., 
H3)) but the method we use suffices for our purposes. 
One of the interesting properties of the formula (|2.8|l is 
that it relates a geometric property of (the n — > oo limit 
of) a bond-diluted graph with the derivative of the Potts 
model, evaluated at a certain temperature as q — > 1, on 
a graph with no bond dilution. Specifically, the mapping 
(|2.ti|) . in conjunction with eq. Q2.8p . formally associates 
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the interval < p < 1 with the interval < v < oo, 
which is the physical range of values of the temperature 
variable v for the ferromagnetic Potts model. In passing, 
we recall that the connection of percolation to the q = 1 
Potts model implies that the percolation transition on 
2D lattices is second-order, in the same universality class 
as the latter model, with the exactly known thermal and 
magnetic critical exponents yt — 3/4 and yt = 91/48 
and thus a = a' = -2/3, (3 = 5/36, 7 = 7' = 43/18, 
5 = 91/5, etc., the values of which were later illuminated 
further by the application of conformal algebra. 

III. SOME BASIC PROPERTIES OF (k) 

As background for our exact results, it will be useful 
to note some basic properties of (k). These can all be 
understood using elementary methods. First, in the limit 
as p — > 0, the only nonzero contribution to the sum in 
the numerator, and to the sum in the denominator, of 
eq. (12.11) arises from the spanning subgraph G' with no 
edges, e(G') = 0, which has n connected components, 
so lim p ^o(^)n = hnip^o(fc) = 1- Second, in the limit 
as p — > 1, the only nonzero contribution to the sum in 
the numerator, and to the sum in the denominator, of eq. 
(|2.1() arises from the term in which the spanning subgraph 
G' is G itself, with k(G) = 1, so lim p ^i(fc)„ = 1/n and 
hence lim p ^i(fc) = 0. Third, for a given {G}, the average 
per-site cluster number (fc) on {G} is a monotonically 
decreasing function of p for p £ [0, 1]. To see this, we 
start with (k) n on a finite G and observe that increasing p 
in the interval [0, 1] decreases the number of components 
k(G). Dividing by n and taking the limit n — > 00 yields 
the corresponding monotonicity result for (k) on {G}. 

Further, for the infinite-length, finite- width strips con- 
sidered here, one can show that (k) is a (real) analytic 
function of p for p £ [0, 1) as follows. We use the general 
relation eq. (|2.8() and observe that the mapping (|2.6() 
formally associates the value p = in the bond perco- 
lation problem with the value v — 0, i.e., infinite tem- 
perature, in the Potts model free energy. Since the free 
energy of a spin model is always analytic near (3 = 0, 
i.e., v = 0, it follows that for any {G}, (k) is analytic 
in the neighborhood of p = 0. As p increases from 
to 1, the value of v corresponding to it via the mapping 
(|2.6() increases from v = to v = 00 (i.e., T = 0). Since 
infinite-length, finite- width lattice strip graphs are quasi- 
one-dimensional, we can use the property that the free 
energy of a quasi-one-dimensional spin system with finite- 
range spin-spin interactions, such as the Potts model, is 
analytic for all finite temperatures. It follows that (k) 
is analytic for < p < 1. In general, for quasi-lD sys- 
tems, the critical percolation probability is p c = 1, so 



that many quantities in ID percolation are nonanalytic 
at p = 1; examples include the average cluster size func- 
tion S(p) = (1 + p)/(l — p), which diverges as p — > 1~, 
and the percolation probability P(p), which is zero for 
< p < 1 and jumps discontinuously to 1 at p = 1 
|2j, |llj, |l7| . This is analogous to the fact that the critical 
temperature for ferromagnetic spin models with short- 
range spin-spin interactions in ID is T c = (v = 00). 
However, for all of the cases that we have calculated, the 
mean cluster number (fc), is analytic at p — 1 as well 
as for p £ [0, 1). This is somewhat similar to the situa- 
tion with the internal energy of a ferromagnetic ID spin 
system with short-range spin-spin interactions, which is 
analytic at T = even though the susceptibility and cor- 
relation lengthy diverge as T — > and the magnetization 
jumps from to 1 at T = 0. This analyticity of (k) for 
p £ [0, 1] for infinite-length, finite-width lattice strips is, 
of course, different from the situation on lattices A of di- 
mensionality d > 2, where (k) generically has a (finite) 
nonanalyticity at p = p c .\- The fact that (k) for the 
infinite-length, finite-width lattice strips is analytic for 
p S [0, 1] does not prevent this quantity from having sin- 
gularities at unphysical values of p, including real values 
outside the interval < p < 1 and complex values. These 
will be discussed below. 

The property that (k) (calculated in the limit as L x — > 
00) is independent of the longitudinal boundary condi- 
tions imposed on the lattice strip follows from the same 
property for the Potts model free energy. As noted above, 
we further expect that, for a given transverse boundary 
condition and for a given p £ [0, 1], the value of (k) for an 
infinite-length strip of the lattice A approaches the cor- 
responding value of (fc) for the (infinite) two-dimensional 
lattice A as L y — > 00. This means that for moderate 
finite widths L y , the resultant (k) can serve as exactly 
calculated approximation to the value of (k) for the cor- 
responding two-dimensional lattice. Of course, this ap- 
proximation refers to the value only, not the full analytic 
structure; as we have shown, for any finite L y , regard- 
less of how large, (k) does not contain the (finite) non- 
analyticity at p c ,A that is present in (k) for the infinite 
two-dimensional lattice. 

IV. RELEVANT RESULTS FOR BOND 
PERCOLATION ON TWO-DIMENSIONAL 
LATTICES 

A. Values of p c ,A 

Here we briefly recall some results that we shall use 
for bond percolation on two-dimensional lattices. For 
a planar lattice G, the (planar) dual, G* is the lattice 
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whose vertices and faces are given, respectively, by the 
faces and vertices of G and whose edges link the faces 
of G. For the thermodynamic limit of a regular two- 
dimensional lattice graph G, the critical probability p c 
for bond percolation satisfies Q p c ,{G} +Pc.{g*} = !• 
In particular, it follows that, since the square (sq) lat- 
tice is self-dual, p c , sq = 1/2. The values of p c are also 
known exactly for the triangular (tri) and honeycomb 
(he) lattices 3]: p c ^ trl = 2sin(7r/18) = 0.347296.. and 
Pchc = 1 - Pc,tn = 0.6527036... 

It is useful to recall that for a general lattice A, p C: \ can 
be obtained from a knowledge of the critical temperature 
of the q — > 1 Potts ferromagnet. The relation l|2.6|) is 
equivalently written as p — v/(l + v). Hence, the critical 
percolation probability is determined as p Cj A = u C) a/(1 + 
v c ,a)> where v Ct \ is the value of v corresponding to the 
phase transition temperature of the Potts ferromagnet 
on A in the limit q — > 1. The value of i> c ,a is known 
exactly for the square, triangular and honeycomb lattices 
(lif . Thus, v c for the square lattice is given by the self- 
duality relation v = Jq; setting q = 1 yields v CtSq = 1 
and substituting this into eq. I|2.6|l yields p c ,sq = 1/2. 
The value of v c jri for the triangular lattice is determined 
as the (unique) real positive solution of the equation Q 
v 3 + 3v 2 — q = 0. Solving this for q — 1 and subsituting 
into eq. (|2.6|) yields the expression above for p c .tri- From 
the duality of the triangular and honeycomb lattices, one 
obtains the critical equation for the latter from the former 
by the inversion map v/^fq —> sfq/v, which yields v 3 — 
3qv — q 2 — 0. Again setting q = 1, calculating the unique 
positive root of this equation, and substituting into <|2.6I) 
gives the expression for p c j lc . 

We next comment on p c for the kagome lattice and 
first recall the definition of this lattice. An Archimedean 
lattice is a uniform tiling of the plane by one or more 
regular polygons with the property that every vertex is 
equivalent to every other vertex. An Archimedean lat- 
tice is specified by the ordered product Y[ vt > where, as 
one makes a small circuit (say clockwise) around any ver- 
tex, one traverses the sequence of polygons pi, with the 
possibility that a given pi occurs ai times in a row (for 
notation, see, e.g., [is))- In this notation, the kagome 
is defined as (3 ■ 6 • 3 • 6). Although the critical bond 
percolation probability p c ,kag is not known exactly, it 
has been determined numerically to high precision |2l| 
as p c ,kag = 0.5244053(3), where the number in parenthe- 
ses is the estimated uncertainty in the last digit. 



expansions around p = and p — 1 . Using results from 
HIH, Ref. [i^] obtained explicit exact values of (k) c for 
bond percolation on three two-dimensional lattices: 



(k) c.sq — 



3 3/2 



and 



(k)c,tri ^ 



= 0.0980762... 



= 0.111844.. 



Petri 



and, using also Q, 



(k)c,hc = 2«fc>c,t« +P 3 c,tri) = 0.07686667. 



(4.1) 



(4.2) 



(4.3) 



(Here we take into account that (k) is defined per site, 
while the quantity nf~ HC in [22| is defined per unit cell 
and there are two sites per unit cell on the honeycomb 
lattice.) Since p c .kag is not known exactly, neither is 
(k) c ,kag- We are not aware of numerical simulations of 
bond percolation on kagome lattices that have yielded an 
approximate value of {k) c ^ a g, or of reasonably long se- 
ries expansions that could be used to obtain an accurate 
estimate of this quantity. 



C. Series Expansions 

Taylor series expansions played an important role in 
the early investigation of critical exponents and tests of 
scaling for percolation. These expansions were performed 
about the points p = (low-density) and p — > 1 (high- 
density), or equivalently, r — > 0, where 



1-p 



(4.4) 



These series expansions have typically been calculated 
for one of two related quantities, the average number, per 
site, of bond clusters, (kbond), or the average number, per 
bond (pb), of bond clusters, (ki, on d)pb- For an (infinite) 
lattice A with coordination number k\, there are ka/2 
bonds per site, so (k bond ) = (n A /2)(k bond ) pb . The aver- 
age number of components (= site clusters) per site, (k), 
differs from the average number of bond clusters per site, 
(kbond) in that (k) counts all components, including iso- 
lated sites, whereas (k bonc i) counts the bond clusters that 
have a nonzero number of bonds and therefore excludes 
isolated sites. Hence, for an (infinite) lattice A, 



B. Values of (fc) c ,A 

Because (k) has only a finite nonanalyticity aip — p c ,A, 
one can obtain rough values of (fc) c ,A from Taylor series 



(fc) A = (kbond) A + (1 ~pY 



(4.5) 



Clearly, (k bo nd) —> rather than 1 as p — > 0, since no 
bond clusters with a nonzero number of bonds are present 
on G in this limit. 
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Early calculations of series for the average bond cluster 
number (kbond) were carried out in Refs. [jj Q and, to 
higher order in 9] for the square, triangular, and honey- 
comb lattices. These have the form 

oo 

(hand) =5> s £> s (r) (4.6) 

8=1 

where r is given by eq. Q4.4jl and the D s are the perimeter 
polynomials. Using the relation (|4. f>|> . it is straighforward 
to calculate series for (k) from those for (kbond)- For our 
present purposes, we shall need only the first few terms 
of these series. For p — > these are 

(k) sq = l-2p + p 4 + 2p 6 - 2p 7 + 7p 8 + 0(p 9 ) (4.7) 
(k) M = l-3p+2p 3 + 3p i + 3p 5 + 3p 6 + 6 P 7 + 0(p 9 ) (4.8) 

(k)hc = 1 - \p + \p 6 + \p W + 0( P n ) . (4.9) 

Thus, for small p, (k) is linear, with slope — K\/2, so 
that (k) decreases more rapidly on a lattice with larger 
coordination number. 

For the p — > 1 series, we have 

(k) sg = r 4 + 2r 6 - 2r 7 + 7r 8 + 0{r 9 ) (4.10) 
(k)tri = r 6 + 3r 10 - 3r n + 2r 12 + 0(r 14 ) (4.11) 

(k) hc = r 3 + |r 4 + V + V + 0(r 7 ) . (4.12) 

These have the general form (k)\ = r KA + ... as p — ► 1, so 
that the larger the coordination number is, the smaller 
(k) in these expansions. 

Because of the significant finite-width effects on the 
lattice strips considered here, one does not expect the 
Taylor series expansions of the exact expressions for (k) 
on these strips to match many orders of the expansions 
for the corresponding two-dimensional lattices. However, 
our exact results can give new insight into one property 
of these series. In early work it was found that the radii 
of convergence of Taylor series expansions around both 
p = and p = 1 were typically set by unphysical singu- 
larities, and these radii of convergence were less than the 
distance from the expansion point to the physical singu- 
larity, p c ,A, for the small-p expansions and r c \ — 1 —p c ,A 
for small-r expansions |9j . (This was also the case for 
three-dimensional lattices [lOj . ) This is reminiscent of 
the situation for low-temperature Taylor series expan- 
sions for various discrete spin models such as the Ising 
model, whose radii of convergence in the expansion vari- 
able z = e~" J were often smaller than the distance from 



the origin to the actual critical point z c , again as a con- 
sequence of complex-temperature singularities j^, E3- 
Although spin models with finite-range spin-spin inter- 
actions only have possible phase transitions at T = 0, 
the complex-temperature singularities of the free energy 
for such quasi-lD systems does exhibit some similari- 
ties with those of 2D systems These complex- 
temperature singularities are associated with complex- 
temperature phase boundaries B, which are the contin- 
uous accumulation set of the (Fisher 0) zeros of the 
partition function. In the same way, we can use our exact 
calculations of (k) for infinite-length, finite-width lattice 
strips to gain some insight into the unphysical singular- 
ities encountered in the above-mentioned series studies 
for percolation. 

V. STRIPS OF THE SQUARE LATTICE 

A. Ly = 1 

The well-known result 

(k) 1D = l-p (5.1) 

for the infinite line can be derived directly using proba- 
bility methods. Here we illustrate how it can be derived 
via eq. (|2.8I) . An elementary calculation yields the Potts 
free energy /(If, q, v) = \n(q + v). Using eq. (|2.8|l yields 
the above result for (fc)i_D. This has the value 1/2 at 
P = Pc,sq (see Table [J. Note that the bond cluster num- 
ber per site in ID is (kbond) — p(l — p), in accord with 
the relation l|4.5|) . 

B. Free Transverse Boundary Conditions 

For L y > 2, we label an infinite- length strip of width 
L y of the lattice A with given transverse boundary con- 
ditions BC y as A, (Ly)BCy- ln particular, the L y = 2 
square-lattice strips with free (F) and periodic (P) trans- 
verse boundary conditions are denoted sq, 2p and sq, 2p. 

1. 2 F 

The free energy of the Potts model for the sq, 2 p strip 
is 

f (sq, 2 F ,q,v) = ~ hi \ sqi 2F,i (5.2) 

where 

^sq.2F,j — ^ (Ts2F ± V R S 2F ) (5.3) 
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with j = 1,2 corresponding to ± and 

T s2F = v 3 + Av 2 + 3qv + q 2 



2F 



Av a 



2„,3 



2 9 u 4 - 2(fv 



I2v 4 



+ 16gw 3 + 13q 2 v 2 + 6q*v + q 



(5.4) 



(5.5) 



In eq. 15.3( 1 only the j = 1 term is relevant for the free 
energy, while the j — 2 term will be discussed below. 

From eq. 1 (5.2( 1 we calculate the average cluster number 
per site 



(k) 



sq,2 F 



(l-p) 2 (2+p-2p 2 
2(1 -p 2 + p 3 ) 



(5.6) 



This is plotted in Fig. ^ together with cluster numbers 
calculated for other strips. At p = p c ,sq = 1/2, this 
average cluster number has the value {k) sq ,2 F = 2/7 ~ 
0.28571. 

From the exact expression for {k) sq ^ F we compute the 
respective Taylor series expansions 

(k) sq , 2F =l-\p+\p i + \p & + 0{p 1 ) for 0(5.7) 

and, in terms of the variable r in eq. I|4.4|l , 

(k) sq ,2 F = ^r 2 + 2r 3 -^r 5 + 0(r 6 ) for r -» . (5.8) 

Thus for this strip (k) is linear for small p and vanishes 
quadratically as p — > 1. As expected for such a small 
width, these series differ from the series for the square 
lattice, although the linear behavior for small p is com- 
mon to both. 

The expression for (k) for this strip has singularities, 
which are simple poles, where the denominator 1 — p 2 + 
p 3 = 0, at 



P ' 



-0.7549, 0.8774 ±0.7449i 



(5.9) 



The first of these poles is the closest to the origin and 
determines the radius of convergence of the small-p Tay- 
lor series in eq. (|5.7|l to be approximately 0.7549. The 
complex pair are the same distance from the point p = 1 
and imply that this series converges for |1 — p\ < 0.7549. 
Thus, although {k) sq ^ F is an analytic function of p for 
p £ [0, 1], the Taylor series expansions about p = and 
p = 1 have radii of convergence less than unity because of 
singularities of this function at real and complex values 
outside the physical interval < p < 1. It is interesting 
that although X sq ,2 F is an algebraic function of v (hence 
p, via eq. (|2.6|l . the resultant expression for {k) sq ^ F is a 
rational function of p. However, this is a consequence of 
the small value of L y . The same comment applies to the 



property that (k) sqt 2 F is meromorphic, i.e., its only sin- 
gularities are simple poles. As will be seen, these features 
are also true of the cluster number (k) for other L y — 2 
strips considered here. 

In this very simple context of a quasi-lD strip, one 
hence gains some insight into the similar influence of un- 
physical singularities in series expansions about p = 
and p = 1 for percolation on higher-dimensional lat- 
tices. To understand these poles more deeply, we ob- 
serve that although the free energy f(sq,2p,q,v) only 
depends on A sg) 2F,i> the partition function for free longi- 
tudinal boundary conditions 45] and q ^ 1 in general is 
a symmetric sum of L^th powers of both of the A M 9F.j's 
for j — 1 and j — 2 (given as eq. (5.17) of Ref. |27|). 
We are interested in the limit q — *■ 1. It is necessary to 
take account of a subtlety concerning the dependence of 
the complex-u phase boundary B of the Potts model, as 
a function of q. In previous work (see eqs. (2.8)-(2.12) 
of Ref. 0, eq. (1.9) of Ref. 0)) we have pointed 
out the noncommutativity at certain special values of q s , 
including q s = 0, 1, namely 



lim lim Z(G, q, v) 

n — >oo q — >q s 



\ 1/n j£ lim lim Z(G,q,v) 1/n 

q—->q s n—*oo 

(5.10) 

and we have noted that because of this noncommuta- 
tivity, for the special set of points q = q s one must 
distinguish between (i) (B({G},q s )) nq , the continuous 
accumulation set of the zeros of Z[G, q, v) obtained by 
first setting q ~ q s and then taking n — > oo, and (ii) 
(£>({G}, q s )) q n, the continuous accumulation set of the 
zeros of Z(G,q,v) obtained by first taking n — > oo, and 
then taking q —t q s . For these special points, 



(B({G},q s )) nq ^(B({G},q s )) 



qn 



(5.11) 



A previous case of this was the q — 2 (Ising) special 
case of the Potts model. Indeed, in that case it was 
noted that B qn does not have the inversion symmetry 
e K — > e~ K that characterizes the Ising model and its 
complex-temperature phase boundar y B na for a bipartite 
lattice (see pp. 396, 433-435 of Ref. |2^V This noncom- 
mutativity is also present at the value q — 1 relevant for 
percolation. If one uses the definition B nq with q = 1 
for percolation, as one uses the definition B nq with q = 2 
for the Ising model, then while B nq is nontrivial for the 
Ising model, B nq is trivial for the percolation problem. 
The reason for this is that if one sets q — 1 first, then, 
from the Hamiltonian definition of the Potts model, since 
the spins are the same on all sites, the spin-spin interac- 
tions on each bond contribute a factor e K to the partition 
function, so one has the elementary result 



Z(G,l,v) = e 



K e(G) 



(v + l) 



e(G) 



(5.12) 
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Substituting v — v p as in eq. i|2.6|) gives 

Z(G,l,v p ) = (l-p)-< G ^ (5.13) 

Evidently, Z(G,l,v p ) has no zeros, so that B nq = in 
the complex p plane. Equivalently, Z(G, 1, v) has only 
a single zero at the point v = — 1, which maps, via eq. 
(|2.6|) . to the circle at infinity in the complex p plane. 
We have noted above that in general, for 5 ^ 1, the 
partition function of the Potts model consists of a sym- 
metric sum of L x 'ih powers of X sq ,2 F ,i and X S q,2 F ,2] if 
one sets q = 1, the coefficient of (\ S q.2 F ,2) Lx vanishes, 
and the (A sgi 2 F ,i) Lx term, with its coefficient, reduces to 
the form (|5.13|) (where in the labelling convention of Ref. 
[27). L x + 1 denotes the number of squares on the sq, 2p 
strip). 

However, the fact that eq. (|2.8|) involves a derivative 
means that it is sensitive to properties of the Potts model 
in the neighborhood of the point p — 1 as well as at this 
point. This suggests that one consider the possible role 
of the locus B qn , although one must use caution in doing 
this because of the noncommutativity discussed above. 
Below we shall use the notation B qn to mean specifically 
the boundary defined for n — > 00 and q — > 1, relevant to 
the percolation problem. 

We find some intriguing connections between the locus 
Bq n and complex-p singularities in (k). Let us first cal- 
culate B qn for the sq, 2 p strip as q — > 1 but with q =/= 1 . 
Evaluating eq. I|5.3|) for q — > 1, we obtain 

A sg ,2F,i = 1 35 ( 5 - 14 ) 



(l-p)< 



and 



A sg ,2F,2 = (1 ^ p)2 ■ (5-15) 

The locus is the set of solutions of the equation of 
degeneracy in magnitude of dominant A's. This locus 
can be seen as a special case of the more general phase 
boundary for the Potts model in the v plane for a fixed q, 
or in the q plane for fixed v. For the present case, since 
there are only two A sg .2Fj's, for j — 1, 2, this equation is 

|A S g,2F,l| = |A s <j,2F,2|j i- e -, 



b 2 (i- P )l = i 



(5.16) 



In terms of the polar coordinates p = pe 10 this equation 



reads p (1 



2 p cos, 9) = 1. The solution is a closed wnere 



egg-shaped curve, shown in Fig. |2J that crosses the real 
p axis at p ~ —0.7549 and p ~ 1.466 and the imaginary 
p axis at p ~ ±0.8688i. This thus constitutes the phase 
boundary in the complex p plane, separating this plane 
into two regions. As follows from the general discussion 
above, the physical interval < p < 1 lies entirely in one 
phase. The three poles of (k) sq _2 F listed in eq. (|5.9() lie 
on this boundary B qn - 



2. 3f,4 f ,5f 

The Potts model free energy / for the infinite-length 
sq, 3f strip was calculated in Ref. [2^]. The free energy 
is given by f(sq,3p,q,v) = (1/3) In X sq .3 F , where X sq .3 F 
is the (maximal) root of an algebraic equation of degree 
4. Because of the complicated nature of the expression 
for this quartic root, we do not present it here. We have 
calculated f(sq, 3f, q, v), and hence (k) sq ^ F , to high pre- 
cision by numerically solving for f(sq, 3f, q, v) for a range 
of values of q near unity, for each value of p, and carrying 
out the differentiation in eq. JUSJ. Although this is nu- 
merical, the computational steps can be carried out with 
almost arbitrarily high precision, so that, in practice, it 
is essentially equivalent to evaluating an explicit exact 
analytic expression. We also apply this procedure for 
larger strip widths, using the exact calculation of / for 
sq,Ap and sq, 5f in Ref- (see also The result- 

ing values of (k) are plotted as functions of p in Fig. ^ 
and the values of (k) at p — p c , S q = 1/2 are listed in Ta- 
ble HJ One could carry out similar calculations of (k) for 
larger values of L y , but our results are sufficient to show 
the nature of the approach of (k) on these infinite-length, 
finite-width strips to the average cluster number for the 
corresponding infinite two-dimensional lattice. Indeed, 
one of the most interesting pieces of information that we 
get from our results - the exact determination of singu- 
larities of (k) in the complex p plane and their effect on 
the radii of convergence of series expansions, can only be 
obtained for strip widths that are small enough so that 
we can get exact explicit analytical forms for (fc). 

C. Periodic Transverse Boundary Conditions 

1. 2 P 

By using periodic transverse boundary conditions, one 
minimizes finite-width effects in this transverse direction. 
We consider first the sq, 2p strip. Note that this strip has 
double transverse bonds. The free energy was computed 
in Ref. [2^] and is given by 

f(sq, 2 P , q, v) = ^ In X sqi2 p,i (5.17) 



A 



s2P 



(5.18) 



with j = 1,2 corresponding to ± and 

T s2P = 6v 2 + 4qv + q 2 + 4w 3 + qv 2 + u 4 (5.19) 

and 

R s2P = ( w 4 + 6v 3 + 8v 2 + 3qv 2 + 6qv + q 2 ) 
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FIG. 1: Plots of (k) (vertical axis) as a function of p £ [0,1] 
(horizontal axis) for infinite-length, finite- width strips of the square 
lattice. The dashed and solid curves refer to free and periodic 
transverse boundary conditions, respectively. For a given p, the 
dashed curves are, in order of descending value of (fc), for lp < 
(L v )p < 5p, and the solid curves are, in the same order, for 2p < 
(L y )p < 5p. 




FIG. 2: Plot of the boundary B qn in the complex p plane for the 
infinite-length sq, 2p lattice strip. Horizontal and vertical axes are 
Re(p) and Im(p). 



The expression (|5.21|) has poles where 1 + p 
ishes, at 



Pi, 2 



i(l±V5)~l 



618 , 



-0.6180 



and where 1 — p + p 2 vanishes, at 



P3A 



1(1 ± V3 



0.5± 0.866i 



p 2 van- 



(5.22) 



(5.23) 



The second and first of these poles are closest to the 
points p = and p = 1 and determine the radii of con- 
vergence of the respective Taylor series expansions about 
these points both to be 0.618. These expansions are 



(k) sq ,2 P = l-2p+^p 2 



2p 4 + 0(p 5 ) 



and 



(k) 



sq,2 r 



r 



2r 4 - 2r 5 + 0(r 6 ) 



(5.24) 



(5.25) 



Note that, in accord with the fact that the coordination 
number of this and any infinite-length lattice strip of the 
square lattice with periodic transverse boundary condi- 
tions is 4, the coefficient of the linear term in the small-p 
expansion is equal to that of the expansion for the full 
square lattice. 

We next discuss the connection of the poles in {k) sq 2 P 
with the locus B qn . Although f(sq, 2p, q, v) only depends 
on the quantity X sq ,2P,i, the Potts model partition func- 
tion for q =/= 1 involves a symmetric sum of i^'th powers 
of both X sq ,2P,i and \ sq ,2P,2 35]. Evaluating \ sq ,2P,j f° r 
q — > 1, we have 



A 



1 



sq,2P,l 



A 



sq,2P,2 



(i-pY 



[i-pf 



(5.26) 



(5.27) 



The locus B qn is the set of solutions of the equation 

|A sq ,2P,l| = |A S( j,2P,2|, i-C 



Ip(l-p)l = l. 



(5.28) 



x (v 4 + 2v 3 + Av 2 - qv 2 + 2qv + q 2 ) .(5.20) 
From this, using eq. I|2.8|l . we calculate 
(l-p) 2 (2-3p 2 + 2p 3 ) 



{k}sq,2p = 



2(1 +p-p 2 )(l -p + p 2 ) 



(5.21) 



At the value of p c , S q = 1/2 for the infinite square lattice 
this has the value (k) sq ^ P = 1/5. 



In terms of the polar coordinates defined above, this 
equation reads p 2 (l + p 2 — 2 p cos 9) = 1. The solution 
forms a closed oval curve in the complex p, shown in Fig. 
13 that crosses the real axis at the points p\,2 in eq. (|5.22|) 
and the imaginary axis at the points p ~ ±0.78615i. As 
in the case of the sq, 2p strip, this curve separates the p 
plane into two regions. All of the four poles of {k) sq ^ P 
given in eqs. I|5.22|l and l|5.23|l lie on this curve B qn . 
This property - that the singularities of (fc) lie on B qn - 
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is analogous to the property that the singularities of ther- 
modynamic functions of spin models lie on the complex- 
temperature phase boundaries for these models, as we 
have studied in earlier work Having pointed 

out the connection between these singularities and the 
locus B qn , we shall, for the strips considered below, just 
summarize the singularities of (k). 




FIG. 3: Plot of the boundary B qn in the complex p plane for the 
infinite-length sq, 2p lattice strip. Horizontal and vertical axes are 
Re(p) and Im(p). 



2. 3p,4 p ,5p 

For the sq, 3p strip, f(sq,3p,q,v) = (1/3) lnA sg ,3 p , 
where A S9) 3 P is the (maximal) root of a cubic equation. 
Although it is possible to display an analytic result for 
(k}sq,3 P j it is sufficiently cumbersome that we do not give 
it here. It is an algebraic, rather than rational, function 
of p. We do display the small-p expansion, which is 



(k)sq. 



1 



3P 



1 - 2p + -p A + p 4 + 0( P b ) 



(5.29) 



was calculated for the sq, 4p and sq, 5p 
, and \ S q,A P and \§ P are roots of equa- 



The free energy 
strips in Ref. [3J 
tions of too high a degree to allow an explicit analytic 
solution. Accordingly, we compute (k) by the numerical 
procedure discussed above. Results are given in Fig. Q] 
and Table H| 

D. Self-Dual Strips of the Square Lattice 

It is of interest to calculate (k) for strips of the square 
lattice that maintain a property of the infinite square 



lattice, namely self-duality. The strips with free and pe- 
riodic transverse boundary conditions considered above 
are not self-dual. However, one can construct a cyclic 
strip that is self-dual by adding a single external site to 
a cyclic square-lattice strip of width L y and then adding 
bonds connecting all of the sites on one side of the strip 
to this single external site. We denote a self-dual (sd) 
strip of this type as sq, (L y ) s d- Before presenting our 
calculations, a remark is in order concerning p c for these 
strips. The physical meaning of p c for a usual infinite 
lattice is, as mentioned before, that for p > p c there ex- 
ists an percolation cluster linking two points that are an 
arbitrarily large distance apart. Now consider the sim- 
plest of the cyclic self-dual lattice graphs, with L y = 1; 
this is a wheel graph, having a rim forming a circuit and 
a central site (~ axle) connected to the sites on the rim 
by L x bonds forming spokes. Evidently, even in the limit 
L x — ► oo, the maximum distance between any two sites 
on this lattice graph is 2 bonds; to get from any site on 
the rim to any other site, one takes a minimum-distance 
route that goes inward along one spoke to the central 
site and out again on another spoke to the other site. 
Similarly, for any finite L y , even as L x — > oo there is 
a maximal finite distance 2L y bonds between any two 
sites. Therefore, although this family of cyclic lattice 
strips does maintain the property of self-duality of the 
infinite square lattice, the notion of a critical p c beyond 
which there is a percolation cluster linking two sites ar- 
bitrarily apart is not applicable to it since no sites are 
arbitrarily far apart. 



I. 1„ 



The free energy is |31|, |3J 

f(sq, l sd ,q,v) = ln\ sd i 

where 

1 



with 



T s di = 3w + q + v 



R sdl = 5v 2 + 2vq + 2v i + q- 2v 2 q + v 4 
From this we calculate 



(5.30) 

(5.31) 

(5.32) 
(5.33) 

(5.34) 



We have (k) sq ^ Bd = 1/6 at p = p c , S q- The mean cluster 
number (k) in eq. I|5.34(l has the following Taylor series 
expansions for p — > and p — > 1 : 

(k)sq,i sd = l-2p+p 3 +p i -p 6 -p 7 +p 9 + 0(p 10 ) (5.35) 
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+ r 4 - r 6 - r 7 + r 9 + 0{r w ) . (5.36) 



VI. STRIPS OF THE TRIANGULAR LATTICE 



One sees that the coefficient of the linear term in the 
small-p expansion correctly matches that of the series for 
the infinite square lattice and the power of the leading- 
order term in p — > 1 expansion is 3, which, although 
not equal to the power 4 in the corresponding expansion 
in eq. ^4.10(1 . is at least closer than the power of 2 for 
the L y = 2 square-lattice strips with free or periodic 
boundary conditions. The poles in eq. H5.34f) at p = 
(1/2) (1 ± v3£) set the radii of convergence of the small-p 
and small-r expansions as unity in both cases, i.e., the 
full physical interval < p < 1. 



A. Free Transverse Boundary Conditions 

1. 2 F 

The free energy for the Potts model on this strip is |2£ 



f(tri,2 F ,q,v) = ^\n\ t2 F 



where 



\t2F 



T t 2F + (3u + v 2 + q)y/Rt2F 



(6.1) 



(6.2) 



with 



t2F 



v 4 + 4iT + 7v 2 + Aqv + q 1 



5 , 4 S 



For these strips, the free energy has the form |31l |3 
f(sq,(L y ) sd ,q,v) = (1/L y ) In \ sq ,(L v ) sd , where \ sq ,(L y ) sd 
are maximal roots of algebraic equationa of degree 5 or 
higher. Hence, it is thus not possible to obtain a closed- 
form analytic solution for this root. We thus follow 
the same high-precision numerical procedure as described 
above. 




and 



Rt2F = q 2 + 2qv - 2qv 2 + 5v 2 + 2v 3 + 



From this we calculate 



(k)tri,2 F 



(I-P) 3 
1 — p + p 2 



(6.3) 



(6.4) 



(6.5) 



Note that this expression for (k) is the same as that for 
the L y = 1 self-dual strip in eq. (|5.34() . This provides an 
illustration of the fact that two different families of lattice 
strips may have the same average cluster number (fc) . We 
plot this cluster number (k) in Fig. together with the 
cluster numbers for the various strips of the triangular 
lattice with greater widths and free or periodic boundary 
conditions. The values of (k) for p = p c ,tr% are listed in 
Table [I] The Taylor series expansions of eq. I|6.5() for 
p — -> and r = 1 — p — > are the same as those of (|5.34H . 



2. 3f,4f,5f 

The free energy / for the strips of the triangular lattice 
with width L y = 3, 4, 5 and free transverse boundary 
conditions were computed in Ref. [3(| (see also [3?|). 
We have used these exact analytic expressions to obtain 
high-precision numerical computations of (k) for these 
strips. 



FIG. 4: Plots of (fc), as a function of p G [0, 1], for infinite-length, 
finite-width self-dual strips of the square lattice. For a given p, in 
order of descending value of (fc), the curves refer to are for l s( j < 
(L y ) s d < 4 S( j. 



B. Periodic Transverse Boundary Conditions 

1. 2 P 

Having explained our calculational method above for 
the square-lattice and previous triangular-lattice strips, 
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0.2 0.4 



0.6 0.8 



FIG. 5: Plots of (fc), as a function of p G [0, 1], for infinite-length, 
finite-width strips of the triangular lattice. The dashed and solid 
curves refer to free and periodic transverse boundary conditions, 
respectively. For a given p, the dashed curves are, in order of 
descending value of (k), for 2p < (L y )p <5f, and the solid curves 
are, in the same order, for 2p < (L y )p < 4p. 



we omit the details for other lattice strips except where 
we have carried out new calculations of Potts model free 
energies. The free energy f(tri,2p,q,v) was calculated 
in Ref. [2^. From it we compute 



{k)tri,2 



(1 - p) 4 (2 + 2p - 7p 2 + 4p 3 - p i + 2p 5 - p 6 ) 



2(1 - 2p 2 



12p 4 + 8p 5 - 2p 6 ) 



(6.6) 



This has the respective Taylor series expansions for p — > 
and p — > 1 : 

(fc)tri,2 P = l-3p+ip 2 + 4p 3 + ^ 4 -10p 5 

- 10/ + O(p 7 ) (6.7) 



{k)tri,2 P = \r A + 2r 6 - 2r 8 + V° + 0(r n ) . (6. 



The poles of (k) occur at 

p ~ -0.3744 , 1.6539 , 

0.1731 ± 0.6306i , 1.1872 ± 0.6924i . (6.9) 

The first two poles, lying on the real p axis, are clos- 
est to the points p — and p = 1 and determine the 
radii of convergence of the series about these points to 
be approximately 0.3744 and 0.6539, respectively. 



Our procedure for strips with greater widths L y > 3 
is as before for the square lattice and (L y )p triangular- 
lattice strips. The small-p series for the tri,3p lattice 
is 



(k) 



tri.3p 



1 - 3p + lp a + 6/ + 0{p 5 ) . (6.10) 



VII. STRIPS OF THE HONEYCOMB LATTICE 

The free energy f(hc,2F,q,v) was calculated in Ref. 
[io| . From it we obtain 

(l-p) 2 (4 + 3p + 2p 2 +p 3 -4p 4 ) 

{k)hc,2 F = 77- 4 . 5A ■ (7.1) 

4(1 — + p a ) 
This has the respective Taylor series expansions 

(k), lc .2 F = 1 - \p + -/ + \ P W - + 0( P U ) (7.2) 



(k)hc 



-r 2 + 3r 3 + 0(r 4 ) 



(7.3) 



The expression (f T . 1 11 has poles at 



p~ -0.8567, -0.15005±0.8975i , 1.0784 ±0.4969i . 

(7.4) 

The first of these is the nearest to the point p = 0, so that 
the small-p Taylor series converges for \p\ ^ 0.8567. The 
last pair of complex-conjugate poles is closest to p = 1, 
so that the series for r — > converges for |1— p\ 0.5031. 

Strips of the honeycomb lattice with other widths and 
boundary conditions are analyzed using the same tech- 
niques as discussed above. The resulting cluster numbers 
(k) are plotted in Fig. |B|and the values for p = p c j LC are 
listed in Table [Q 

VIII. STRIPS OF THE KAGOME LATTICE 



2 F 



For the purpose of obtaining (fc), wc have carried out 
a calculation of the free energy of the Potts model on the 
2p strip of the kagome lattice. The results are sufficiently 
lengthy that we list them in the appendix. From these 
we calculate 



(k)ka 



2F 



S,2f 



D 



k2F 



(8.1) 



where 



N k2F = (1 - p) 2 (5 + 2p - p 2 - 2p 3 - 8p 4 - 16p 5 
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FIG. 6: Plots of (fc), as a function of p G [0, 1], for infinite-length, 
finite-width strips of the honeycomb lattice. The dashed and solid 
curves refer to free and periodic transverse boundary conditions, 
respectively. For a given p, the dashed curves are, in order of 
descending value of (k), for 2p < (L y )p < 5p, and the solid curve 
is for 4p. 



FIG. 7: Plots of (fc), as a function of p G [0, 1], for infinite-length, 
finite-width strips of the kagome lattice. The dashed and solid 
curves refer to free and periodic transverse boundary conditions, 
respectively. For a given p, the dashed curves are, in order of 
descending value of (k), for 2p < (L v )p < 3p, and the solid curve 
is for 2p. 



-43/ - 26p 7 - 2p 8 + 10p 9 - 3p 10 



(8.2) 



and 



D 



k2F 



5(1 - p 4 - 2p 5 + 10p 6 - 10p 7 + 3p 8 ) . (8.3) 



We plot (k) kag: 2 F in Fig. □ At the value p = p cMg , 
(k)kag,2 F nas the approximate value 0.22918. The (fc) for 
this 2 p Kagome strip has the following expansions in the 
vicinity of p = and p = 1 : 



(k) ka 



S,2f 



-P 



+ -P 3 + -=p 6 + 0( P 7 ) (8.4) 



1 
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(k) kag .2 F = ~r + ~r A + — r 4 - 7 r 5 + 0(r b ) . (8.5) 
5 5 5 5 

The cluster number l|8.1|> has poles at 

p ~ -0.5470 ± 0.2862i , -0.0363 ± 0.6583i , 



0.7772 ± 0.5605i , 1.4728 ± 0.1486i 



(8.6) 



Of these, the first and last complex-conjugate pairs are 
closest to p — and p = 1, respectively, and deter- 
mine the radii of convergence of the Taylor series ex- 
pansions about these points to be approximately 0.6174 
and 0.4956. 

We have also calculated (k) for the 3f strip of the 
kagome lattice; this is plotted in Fig. [7\ 



B. 2 P 

For the 2p strip of the kagome lattice we find 

N k2P 



(k) ka 



g^p 



D 



k2P 



(8.7) 



where 

N k2 p = (1 - p) 4 (6 + I2p + Up 2 + 4p 3 - 25p 4 - 108p 5 
+ 16p 6 + 472p 7 - 706p 8 + 320p 9 + 286p 10 - 352p n 
-194p 12 + 360p 13 + 120p 14 - 340p 15 + 65p 16 + 136p 17 
-96p 18 + 24p 19 - 2p 20 ) (8.8) 
and 

D k2P = 6(1 - 2p 4 - 8p 5 + 32p 6 + 40p 7 - 268p 8 
+424p 9 - 320p 10 + 120p n - 18p 12 ) . (8.9) 

A plot is given in Fig. {7\ The cluster number (k) kag- 2 P 
has poles at 

p ~ -0.4660 , 1.6556 , -0.3443 ± 0.2919i , 

-0.0057 ± 0.475H , 0.5325 ± 0.48455i , 

1.0776 ± 0.4384i , 1.4785 ± 0.2140i (8.10) 
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The first complex-conjugate pair is the nearest to the 
origin and sets the radius of convergence of the small-p 
Taylor series expansion of (k)k a g,2 P as 0.4514, while the 
second-to-last complex-conjugate pair is closest to the 
point p = 1 and determines the radius of convergence of 
the series expansion about this point to be 0.4453, to the 
stated accuracy. 

To our knowledge, it is not known what the value of 
(k) is for the (infinite) kagome lattice at the numerically 
determined critical percolation probability p c ,kag- As- 
suming that, for a given set of transverse boundary con- 
ditions and a given p G (0,1), (k) is a monotonically 
decreasing function of the strip width L y for this lattice, 
as we find for other lattice strips, our results yield the 
upper bound (k)kag < (k)kag,2 P — 0.11149 at the value 
V = Pckag given above. Here we use the result for the 2p 
strip since it is lower than the result for the 2p and 3f 
strips. 

IX. DISCUSSION 

We first introduce a notion of effective coordination 
number. For a graph G the degree of a vertex is the 
number of bonds connected to this vertex. A K-regular 
graph is a graph in which all of the vertices have the 
same degree, k. Whether a given lattice strip graph 
is K-regular depends on the longitudinal and transverse 
boundary conditions; for example, it is k regular if one 
uses toroidal (doubly periodic) boundary conditions. In 
the limit L x — > oo, since the longitudinal boundary con- 
ditions do not affect the free energy f({G}, q, v), we need 
only consider the effect of the transverse boundary con- 
ditions. The effective coordination number is 

K eff ({G}) = lim ^ . (9.1) 

n— >oo n(G) 

Clearly K e ff = n for a regular lattice. For regular lat- 
tice strips with periodic transverse boundary conditions, 
the value of n e ff is the same as the value for the cor- 
responding two-dimensional lattice. For strips with free 
transverse boundary conditions, we have 

Keff {A,(L y ) F )= KA (l--^J (9.2) 

where ka = 4, 6, 3 for A = sq, tri, he and 

1 2 1 . . 

a sq = 2 ' a *« = 3 ' a hc = 3 • (9-3) 

For the cyclic self-dual strips of the square lattice, the 
single external vertex connected to each of the sites on 
one side of the strip has a degree L x that diverges as 
L x — > oo. The L x (L y — 1) interior vertices have degree 4, 



TABLE I: Values of (k) on infinite-length strips of the lattices A 
(where sq, tri, he, kag denote square, triangular, honeycomb, and 
kagome) of finite width L y at p = p c Aj The transverse boundary 
conditions (BC' y ) are F and P for free and periodic, respectively. 
We also include results for the self-dual (sd) strips of the square 
lattice. The effective coordination number K e tt is defined in cq. 
19.11 . The L y = oo values of (k) p — P( . A are the values for the two- 
dimensional lattice A, and the dashes for these entries indicate that 
they do not depend on BC y . 
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while the L x vertices on the rim have degree 3. Together, 
these lead, in the limit L x — > oo, to the result K S q,sd = 
4. Finally, for the kagome strips with free transverse 
boundary conditions 

Keff (kag, (L y ) F ) = 4(l - 3L 1 _ J (9.4) 
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while for the kagome strips with periodic transverse 
boundary conditions, K e ff = 4, the same value as for 
the infinite two-dimensional kagome lattice. 

From our calculations we find a number of generic fea- 
tures: 

• We have shown that (k) is a (real) analytic function 
of p in the interval < p < 1. At the critical 
percolation probability p = 1 for these quasi- ID 
strips, our exact results for (k) are also analytic, 
although some other quantities in percolation, such 
as the percolation probability P(p) and the cluster 
size S(j>) are not, as is evident from the well-known 
ID case. 

• As the curves in the figures show, with an increase 
in strip width L v , (k) is consistent with approach- 
ing a limiting function of p. This is in accord with 
one's expectation. 

• For a given p in the interval between and 1, and 
for a given type of lattice strip, as the width L y in- 
creases, (k) decreases, so that the approach to the 
asymptotic value for the 2D lattice is from above, 
in the cases that we have computed. For strips with 
free transverse boundary conditions, increasing L y 
increases K e //, so the decrease of (k) is associated 
with an increase in the effective coordination num- 
ber. This is reasonable, since, heuristically, for a 
fixed value of p, there is a greater probability of 
having a percolating cluster on a lattice of higher 
coordination number, so that more sites are part 
of this cluster and there are fewer separate clus- 
ters per site. This is also reflected in the mono- 
tonic decrease of p c ,h with increasing k\ for most 
higher-dimensional lattices. (However, we recall 
that counterexamples to this general monotonic de- 
crease of (k) with increasing coordination number 
are known |23l l2fij.) 

For strips with periodic transverse boundary condi- 
tions, the decrease of (k) at a fixed p with increasing 
width L y is not associated with an increase in K e //, 
since re e // is constant for these strips (and equal to 
the two-dimensional value) ; here one may interpret 
the decrease as being simply due to a reduction in 
the finite-width effects that enables the percolation 
quantities to approach their two-dimensional val- 
ues. 

• For a given lattice type, we find some examples 
where the curve for (k) calculated on a strip of 
width L y with periodic transverse boundary con- 
ditions will cross the curve for (k) for the same lat- 
tice and a different L v and free transverse bound- 



ary conditions. For example, as is evident in Fig. 
1, the curve for (k) on the sq,2p strip lies below 
those for (k) on the sq, {L v )f strips at small p, but 
sequentially crosses the latter as p increases and lies 
above them (except for L y = 1,2) as p — > l - . Sim- 
ilar behavior is observed, e.g., on the strips of the 
triangular lattice. These also constitute examples 
of how (k) calculated on a strip with a larger value 
of K e f f than that of another strip can be larger than 
(k) for the latter strip. For instance, K e ff = K = 4 
for the sq,2p strip, which is larger than the value 
K e ff = 3.6 for the sq,5p strip; however, (k) on 
the former strip is larger than (k) on the latter for 
p <^ 0.36. This dependence on transverse boundary 
conditions is consistent with disappearing as the 
strip width L y — > oo, consistent with the approach 
to a single limiting function (k) for the correspond- 
ing 2D lattice. Although we have not proved rigor- 
ously that the function (k) obtained via this limit- 
ing sequence (taking L x — > oo first and then taking 
L y — > oo) is identical to the function (k) obtained 
via the usual two-dimensional thermodynamic limit 
(L x — ► oo, Ly — > oo with L y /L x a nonzero finite 
number) , this conclusion is consistent with our find- 
ings. 

• We have used the values of (k) at p = p c ^\ as a 
measure of how rapidly, for a given p, the cluster 
number calculated on infinite-length, finite-width 
strips approaches the value for the two-dimensional 
lattice. These values are listed in Table H] Even for 
the modest strip widths considered here, one sees 
that (i) these values approach the known values of 
(k) on the corresponding two-dimensional lattices 
reasonably quickly, and (ii) this approach is more 
rapid when one uses periodic transverse boundary 
conditions, as is expected, since the latter mini- 
mize finite- width effects. For example, for the strip 
of the square lattice with L y = 5 and periodic 
transverse boundary conditions, (fc) evaluated at 
P = Pc, sq is about 15 % larger than the value 14. 1|) 
for the square lattice, while (k) for the tri,&p and 
hc,4p strips, evaluated at the respective p c ,tri and 
p c ,hc, are both about 17 % larger than the corre- 
sponding values H4.2f> and (|4.3|l for the triangular 
and honeycomb lattices. 

• We find that for these strips, the small-p series ex- 
pansions of (k) have the leading terms 

<*> = 1-(^)p + - ( 9 - 5 ) 

which are analogous to the structure that these se- 
ries have for regular lattices of dimension d > 2. 



Higher-order terms in the series for the strips of 
small widths are not expected to coincide with 
those in the series for the two-dimensional lattices, 
and one sees that they do not. 

• An interesting output of our analysis is the exact 
determination, for various infinite-length, finite- 
width strips, of the singularities of (k) in the com- 
plex p plane. As we have shown, for many strips 
these (real and/or complex) singularities outside 
the physical interval [0, 1] occur sufficiently close 
to the points p — and p = 1 that they render the 
radii of convergence of the respective Taylor series 
expansions about these points less than unity, al- 
though the actual functions (k) themselves are an- 
alytic functions on p 6 [0,1]. Although the strip 
widths are probably too small to justify a detailed 
comparison with unphysical singularities for per- 
colation quantities in two dimensions, this generic 
property - the presence of unphysical singularities 
that determine the radii of the Taylor series expan- 
sions about the points p = and p = 1 to be less 
than p c for the given type of lattice - is similar to 
what was found in analyses of series for theperco- 
lation problem on two-dimensional lattices 9] (and 
three-dimensional lattices [lCj). 

• Finally, we have discussed how, for a given infinite- 
length, finite- width strip, the unphysical singulari- 
ties have a connection with the locus B qn , which is 
the continuous accumulation set of the zeros of the 
Potts model partition function in the p (or equiva- 
lently the v) plane obtained by first letting n — ► oo 
and then q — > 1. In particular, we find that these 
unphysical singularities lie on B qn . The noncom- 
mutativity of eq. (|5.11|l analyzed in Ref. plays 
a crucial role here, since B nq , obtained by first let- 
ting q — > 1 and then n — > oo, is trivial. Our results 
motivate further study on this topic. 

X. CONCLUSIONS 

In summary, we have presented exact calculations of 
the average cluster number per site (k) for the bond per- 
colation problem on infinite-length, finite- width strips of 
the square, triangular, honeycomb, and kagome lattices, 
with both free and periodic transverse boundary condi- 
tions. We believe that these results are a useful exten- 
sion beyond the one-dimensional result toward two di- 
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mensions and provide insight into the form of (k) as a 
function of the bond occupation probability p. 
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XI. APPENDIX 

In this appendix we give the free energy for the Potts 
model on the 2p strip of the kagome lattice. We find 

f(kag,2 Fl q,v) = \\n\ k2F (H- 1 ) 
5 

where 

Xk2F = \[T k2F + ^/Rk2F] (H.2) 

with 

T k2F = v 8 + 8v 7 + v 6 q + 29v 6 + 20 A + lOuV 
+ 2v 3 q 3 + A2v 5 + 61 A + 54«V + 28 A 3 
+ 8vq 4 + q 5 (11.3) 

and 

Rk2F = v 16 + 16v 15 + 2v u q + 114A + 32v 13 q 
-3v 12 q 2 - Av n q 3 + 484t; 13 + 288v 12 q + 52v n q 2 
-2&v 10 q 3 - 12«V ~ 2wV + 1329v 12 + 1572v n q 
+ 1098i; 1 V + 520A 3 + 192A 4 + 48«V + 6A 6 
+2196w n + 4350w 10 g + 5196A 2 + 4344vV 
+2628v 7 q 4 + 1114vV + 312vV + 52v 4 q 7 + 4v 3 q s 
+ 1620w 10 + 4572v 9 g + 7413 A 2 + &284v 7 q 3 + 6732v 6 g 4 

+4028A 5 + 1766 A 6 + 556 A 7 + 120A 8 + 16V 

+<Z 10 . (11.4) 

(It can be checked that the v = — 1 special case of this / 
coincides with the degeneracy per site, W, in Refs. [471].') 
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